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We study the impact of a leptophobic Z' gauge boson on the C\ q and C2 9 parameters that describe 
the low-energy, parity- violating electron-quark neutral current interaction. We complement previous 
work by including the penguin-like vertex corrections, thereby completing the analysis of one-loop 
calculation up to 0{m 2 q , /M%,) terms. We analyze the sensitivity of these probes to the different 
couplings Z'uq (q = u, c, t) and Z'dq (q = d, s, fa), in a model-independent way that can be applied 
to any specific Z' scenario. We show that constraints from neutral kaon and heavy flavor studies 
preclude significant contributions from flavor non-diagonal couplings except for those involving top 
quarks. We apply our results to a light Z' with flavor diagonal couplings to up or down quarks, a 
scenario proposed in the literature to explained the CDF W plus di-jet anomaly. We find that such 
a particle would not affect the C\ q coefficients, but it would have a sizable impact on C*2 9 couplings 
that can be probed by future measurements of parity- violating deep inelastic scattering of polarized 
electrons off of deuterium. 



Despite the impressive success of the Standard Model 
(SM), it is widely believed that it is not the final the- 
ory and that some extension is needed. One of the sim- 
plest modifications is the addition of an extra U(l) gauge 
group, with the consequent addition of a Z' gauge boson. 
This extension is present in many scenarios for physics 
beyond the Standard Model (BSM), and the details of 
the Z' interactions with fermions are completely model- 
dependent. 

The possibility of a leptophobic Z' boson, both with 
diagonal and non-diagonal couplings in flavor space, has 
received some attention recently in the literature as a 
possible explanation for some experimental anomalies [1- 
5]. Interestingly, this kind of Z' scenario is difficult to 
probe with colliders due to the QCD background, and 
the limits on them are not strong. 

In this paper we analyze the possibility of indirectly 
probing these models using parity-violating (PV) deep 
inelastic scattering (DIS) of electrons off of deuterium. 
The measurement of this PV asymmetry has been car- 
ried out recently at Jefferson Lab (JLab) [6], and more 
precise measurements are planned for the JLab 12 GeV 
program [7, 8] and possibly for a future Electron Ion Col- 
lider (EIC) [9]. 

The low-energy effective Lagrangian describing the PV 
interaction of an electron and a light quark q = u, d is 

£pv = —7= ^2 {Ciqin^^iqj^q) + C 2q (e^e)(qj t ,j 5 q)} , 

(1) 

where Gf is the Fermi constant that can be determined 
from the muon lifetime. At tree-level in the SM, these 
couplings are given by 
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C lq = -/f + 2Q g sin 2 
C 2q = If (-1 + 4 sin 2 W ) 



(2) 
(3) 



the weak mixing angle. The 0(a) electroweak radiative 
corrections can be found in Refs. [10-12], and the the- 
oretical uncertainty in C\ q . 2q is below the percent level. 
The coefficients C\ q have been extracted with high pre- 
cision [13] from atomic PV in Cesium [14, 15] and will 
be determined also through PV elastic ep scattering [16]. 
On the other hand, the current experimental determi- 
nations of the C 2q coefficients are much less precise (cf. 
Rcf. [12]), although this situation will change in the near 
future due to PV-DIS studies. We will see in this paper 
how C 2q are sensitive to some BSM models that do not 
affect the C\ q coefficients, showing their complementarity 
as BSM precision probes. 

More specifically, the parity-violating eD asymmetry 
depends on the Ci q coefficients as follows 
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l + (l-y) 2 J' (4) 
(2C lu - C ld ) [1 + Ri] , (5) 

(6) 



- (2C 2u - C 2d ) [1 + R 2 ] , 



where Q q and /| are the electric charge and the third 
component of weak isospin for the quark q, and #w is 



where Q 2 = \q\ 2 — q 2 is the four- momentum transfer 
squared, y is the fractional energy transfer from the elec- 
tron to the hadrons and R\. 2 include different hadronic 
corrections (see Ref. [17] for more details and references). 
The SOLID experiment at JLab will measure this asym- 
metry with a statistical precision better than 0.5% for 
a number of bins with 0.3 < x < 0.75, 4 < Q 2 < 10 
(GeV/c) 2 , and y ~ 1 [8]. SOLID has a significant sen- 
sitivity to the C 2q coefficients due to the large value of 
y, and is projected to determine 2C 2u — C 2 d with an er- 
ror ±0.0083 [18], which could be improved by an EIC 
measurement. 

A leptophobic Z' with vector and axial-vector cou- 
plings to quarks will in general modify the value of the 
Ci q coefficients, and thus it will have measurable conse- 
quences on the PV eD asymmetry. Due to its leptophobic 
nature the modifications can be classified as vertex and 
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FIG. 1: Classes of corrections to the eq scattering due to a 
leptophobic Z': 7 — Z' (or Z — Z') mixing and vertex cor- 
rections in the quark side of the diagram. The external leg 
(quark self-energy) corrections are not shown. 

external leg (quark propagator) corrections in the quark 
side of the diagram and 7 or Z mixing with the new neu- 
tral gauge boson Z', as shown diagrammatically in Fig. 1. 
The latter effect was studied in a recent publication [20] 
that we will complement in this work by analyzing the 
effect of the vertex correction and external leg correction 
(not shown). 

We will follow a phenomenological, bottom-up ap- 
proach, where the Z' mass and its couplings to quarks 
are free parameters. Specifically, we will use the following 
Lagrangian to describe the interaction of the leptophobic 
Z' boson with quarks 

C{x) = g'Z'^ J2 ^ (q'u + Q'iH) 1j • (7) 

ij 

where g' is the new gauge coupling constant and 
are the flavor-dependent vector (axial-vector) quark cou- 
plings, that satisfy Q^^* = Q^ A ^ to ensure the her- 
miticity of the Lagrangian. In this expression the quark 
fields are the mass eigenstates, i.e. the CKM rotation 
has already been performed. 

Both the Z and 7 couplings to quarks receive correc- 
tions from one-loop diagrams containing the new Z' bo- 
son. In turn, these vertex corrections (and the related 
external leg corrections) generate new contributions to 
the PV Ci q coefficients through penguin-like diagrams, 
as shown in Fig. 2. Although we have of course included 
the external leg corrections in our computation, we have 
omitted them from Fig. 2 for the sake of brevity. 
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FIG. 2: New penguin diagram built with the Z' correction to 
the yqq vertex, that generates a non-zero BSM contribution 
to C 2q - 

If the Z' is much heavier than the quark j running in 
the loop, we can neglect terms of order m|/M|, . In this 
case, the sum of vertex and external leg corrections to 



both the 7<7<7 and Zqq interactions vanish in the Q 2 — > 
limit. The leading non- vanishing correction that trans- 
forms as an axial vector coupling occurs at second or- 
der in the four-momentum transfer. For the Zqq in- 
teraction, the resulting contribution to the four-fermion 
(e-f^e) (97^759) operator is suppressed by Q 2 /M%, mak- 
ing it negligible at the kinematics of the future PV-DIS 
experiments. In contrast, the 0(Q 2 ) contribution to the 
axial vector jqq vertex, commonly referred to as the 
"anapole moment" [19], compensates for the 1/Q 2 aris- 
ing from the photon propagator in the penguin graph of 
Fig. 2, leading to the following Q 2 -independent contribu- 
tion to the Ci q : 

SC lq = , (8) 
SC -- 2Q «a ,2 s 2 c 2 M2z O' v O' A ( l W m i) 

where Cw and sw are the cosine and sine of the weak 
mixing angle 9yj. In this expression and in the remainder 

of this article we assume the couplings Q^ A ^ to be real 
for the sake of simplicity. Since the photon does not 
have any axial-vector coupling, the C\ q coefficients do 
not receive any BSM contribution. On the other hand, 
the Z' contribution to the anapole moment of the quarks 
generates the above-given non-zero 5C2 q - 

It is worth stressing that these results are model- 
independent, in the sense that they do not depend on 
the details of the UV completion of the Z' model but 
only on its mass and its couplings to quarks. 

One might be interested in the same expression for the 
case when the mass of the quark running in the loop is not 
negligible compared with the Z' mass, as would happen 
with a 0(100) GcV Z' particle with non-diagonal cou- 
plings to up and top quarks. In this case, however, this 
result depends on the UV completion of the theory and 
hence a modcl-indcpcndcnt evaluation is no longer possi- 
ble. In particular, the Z' contribution to the Z or 7 ver- 
tex correction shows a gauge dependence that will be can- 
celled by other new contributions in the more complete 
theory where the Z' is embedded 1 . This issue was also 
discussed in Ref. [21], where the 0(m 2 /M§,) effect on the 
Ci q coefficient due to a light Z' boson with non-diagonal 
flavor coupling to up and top quarks was studied in a 
specific model where the new neutral gauge bosons are 
introduced through a non-Abelian flavor gauge symmetry 
[22] . Models where the Z' has a fairly large non-diagonal 
coupling between up and top quarks have received some 
attention recently [23] since they could accommodate the 
large measured top quark forward-backward asymmetry 
A l FB at the Tevatron [24] . 



Notice that this situation is the same in the Standard Model if 
we substitute the Z' by the electroweak Z boson. In that case 
the gauge dependence is canceled by one-loop diagrams with the 
scalar degrees of freedom coming from the Higgs sector. 
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It is useful to notice at this point that a Z' boson with 
up-charm, down-strange or down-bottom flavor changing 
couplings would contribute at tree- level to the D° — D°, 
K° — K° or B® — mass splittings respectively. Namely, 
if we assume small Z — Z' mixing we find [25] 
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mp 



g'Fp 
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(9) 



(o;/) 2 (fcp + 2)-(o;j) 2 (fc P -2) 



where mp and Fp are the mass and decay constants of 
the P meson (P — D 7 K,Bd); i,j denote the valence 
quarks in P; and kp is given by 



(mi +m 3 ) 2 



Thus, we see that unless an unnatural cancellation be- 
tween the V and A terms in Eq. (9) occurs, this result 
can be used to set strong bounds on these flavor non- 
diagonal couplings. 

It has been shown experimentally that these three mass 
splittings Ams,D,K are non-zero [12]. These measure- 
ments are actually very precise in the B and K systems, 
where the relative errors are at the 0.1-1.0% level. How- 
ever, the determination of the mass splittings in the Stan- 
dard Model is quite complicated, especially in the D and 
K cases due to long-distance effects, which opens the 
possibility for large NP contribution. Nonetheless, even 
in the extreme scenario where the BSM contribution to 
Amp is of the same size of the experimental value, we 
still find that 



9 '\q7< a ^ 
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< 10~ 3 



(11) 



for the three cases — (u,c),(d,s) 7 (d,b) 2 . Such a 

strong constraint makes their effect on the PV coefficients 
C 2q completely negligible: 



\SC 2q \ « \Q q \ 10" 



6 ~ l0g Ml 



(12) 



We can therefore safely neglect the flavor non-diagonal 
terms in our Ci q expressions, except for the up-top cou- 
pling that we will keep. After some trivial manipulations 



we can then write 
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(14) 



(15) 



(16) 



where we have replaced the light quark masses m u ^d ap- 
pearing in Eq. (9) by a hadronic scale of 1 GeV, as a 
conservative cut-off, since non-perturbativc effects bc- 
(10) come large for smaller energies. The B^ t ' 2 ^ are model- 
dependent quantities with the following asymptotic be- 
havior in the limit M z > mt- 
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(17) 
(18) 



As shown in Fig. 1, a leptophobic Z' also modifies the 
C 2q coefficients through photon-Z' mixing, although it 
cannot affect C\ q as the photon does not have the tree- 
level axial coupling to the electron needed to generate the 
(e7 M 75e) (jti^q) operator. On the other hand the Z — Z' 
mixing diagram does contribute to both C\ q and C 2q co- 
efficients, but the associated mixing angle is constrained 
to be so small [27] that this contribution can be safely 
neglected. 

The 7-Z' correction to C 2q was calculated recently with 
the following result [20] 



5C- 



2q 



9tt 2 



n' 2 <? r 2 
9 s w c w 



Mz 
Mz 



Q j x 



(19) 



10.54 (2Q' u v u - Q' d v d ) - 8.65 Q'/ s + 
12.43 QZ - 4.44 Q'Z - 1.89 qZ 



where we have conveniently chosen a form similar to 
Eqs. (15) and (16) to ease the comparison. 

We observe a complementarity between both contribu- 
tions to C 2q (vertex and propagator corrections): 

• Non-diagonal flavor couplings do not contribute to 
the 7 — Z' (or Z — Z') mixing since the SM neutral 
current couplings are flavor diagonal, but they do 
contribute to the vertex corrections. This is partic- 
ularly relevant for a Z' that couples to up and top 
quarks; 



2 For a more detailed analysis of the effect of a non-universal Z' 
in meson- ant imeson mixing see Ref. [26] and reference therein. 



The dependence on the Z' mass is different, due 
to the presence of a logarithmic term in the vertex 
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correction. However, we wrote Eqs. (15) and (16) 
in such a way that the logarithm contribution is 
subdominant for Z' masses in the 100-1000 GeV 
range; 

Moreover, in the case of flavor diagonal couplings: 

• Obviously the corrections to C 2q are zero unless the 
Z' has axial- vector couplings to the qq pair (Q q A ). 
However, for the vertex correction we need also a 
non-zero vector coupling to the same qq pair (Q q ^), 
whereas in the 7 — Z' mixing the vector coupling 
can be to a different quark pair {Q q Y q ,); 

• The 7 — Z 1 mixing correction is generally larger in 
magnitude than the vertex correction, unless some 
cancellation between the different terms in Eq. (19) 
takes place. 

Therefore, our final result for the Ci q coefficients is 
given by the sum of the vertex corrections of Eqs. (13)- 
(16) and 7 — Z 1 mixing correction given in Eq. (19). In 
particular, if the Z' couplings are flavor universal and 
diagonal these expressions take the following form 
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(20) 



(21) 



-26.92 + log 



M, 



M, 



23.63 Q' d v d 



For the sake of illustration, we now apply our results to 
two specific Z' scenarios that have been studied recently 
in the literature as possible explanations for the W ± +jj 
excess observed by CDF Collaboration in the study of 
events with a lepton, missing transverse energy and a 
pair of hadronic jets [29] 

The first case is given by the leptophobic Eq GUT 
scenario outlined in Ref. [28], where the Z' couplings to 
quarks are flavor diagonal and universal, with Q^ u = 1/6, 
Q'± = 1/2, Q'/ d = -1/3 and = 0. In this model the 
gauge coupling g' and the mass of the Z' boson are not 
fixed, but it is illustrative to use g' <~ 0.6 and Mz' ~ 
150 GeV, since these values can accommodate the CDF 
W ± + jj excess, as shown in Ref. [2] . Hence we find in 
this benchmark scenario at Q 2 — the following result 



6C- 
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150 GeV s 
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M z , J ^0.6 
M z , 
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5C 2d = 0. 



(22) 
(23) 



the conclusions of Ref. [20] are not substantially modi- 
fied by the inclusion of the vertex correction that only 
decreases the result by ~ 7% 3 . It must be emphasized, 
however, that this is not true in general and one could 
have Z' models were the vertex correction is the domi- 
nant. Nevertheless these cases are less interesting phe- 
nomenologically, since the value of SC 2u is then smaller. 
Likewise the absence of a C 2 d correction is due to the 
zero value of the down quarks axial-vector charge {Q d A ) 
in our benchmark Z' scenario and hence is not a model- 
independent feature. 

We note the results shown so far have been obtained 
at Q 2 — 0, whereas the future PVDIS experiments will 
take place at a finite value of Q 2 . In the case of the 
SOLID experiment the Q 2 will be in the range 4-10 
(GeV/c) 2 . In Ref. [20] it was found that this finite-Q 2 
effect reduces the absolute value of SC 2u by <~ 25% (~ 
30%) at the lower (upper) end of the kinematic range, a 
modification that can also be applied also to our results 
for this first example since the final result is dominated 
by the propagator corrections. 

A deviation of this size (2C 2u — C 2c i ~ 0.02) from the 
Standard Model prediction will be accessible to SOLID 
and the EIC, therefore probing this class of Z' models 
and evaluating the presence of axial- vector couplings to 
quarks. 

As a second example, we discussed now the phe- 
nomenological model introduced in Ref. [3] to explain the 
CDF W ± +jj excess [29], through a light leptophobic Z' 
that couples only to first generation quarks. 

As shown in Ref. [3], the collider bounds on such a Z' 
are not strong, and in fact for a light Z' (below 300 GeV) 
the strongest bounds come from the dijet searches at the 
UA2 experiment [30] (pp — >• Z' — >• qq) . At the Tevatron, 
in contrast, the QCD background become too large com- 
pared with the signal. The comparison of these bounds, 
extracted in Ref. [3] , and the bounds that can be derived 
from the measurements of the PV eeqq coefficients C 2q 
are shown in Fig. 3. As we can see, a future determi- 
nation of the quantity 2C 2u — C 2 d with an error in the 
10~ 2 — 10~ 3 range would probe a comparable portion of 
the parameter space as the UA2 measurements. To facil- 
itate the comparison, we follow in this plot the notation 
of Ref. [3] , where the Z' is assume to couple either to left- 
or right-handed quarks, but not both at the same time, 
and where their coupling constant g qq z' is related to our 
coupling constants as follows 
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(24) 



Here, the plus (minus) sign corresponds to the case 
when the Z' boson couples to right-handed (left-handed) 
quarks. 



This correction to the C 2u coefficient, representing ap- 
proximately a 30% correction to the SM value, is com- 
pletely dominated by the Z' — 7 mixing term, and thus 



3 The overall sign of the correction SC2 U is however the opposite 
than in Ref. [20], due to a typo in the sign of the axial- vector 
coupling of up quarks in that work. 
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FIG. 3: Regions of the g qq z' — M z < plane probed by the UA2 
experiment as extracted in Ref. [3], and by future measure- 
ments of the quantity 2C2u — Cid with different levels of pre- 
cision. The upper and lower plots show the result for a Z 1 
that only couples to u and d quarks respectively. 

In particular, it was shown in Ref. [3] that a Z' with 
a mass of 150 GeV and couplings to left-handed first- 
generation quarks could explain the CDF W +jj excess. 
Applying our results, this Z' would produce the following 
correction to 2C 2u — C 2 d 

S(2C 2u -C 2d ) = -0.0l2[4g 2 uuZ ,+g 2 ddz , 

-5.3 g U uZ'9ddZ'\ ■ (25) 



As shown in Fig. 4 of Ref. [3] , the values of the couplings 
9uuZ' and gddz 1 needed to explain the W ± + jj excess 
are not fixed separately. However we can use the point 
(fjuuZ' , 9ddZ') ~ (0.25,0.00) as an illustrative example, 
and in that case we find 8 (2C 2u — C 2 d) ~ —0.003, an ef- 
fect smaller than in the Eq GUT scenario, but still above 
the permil level. Probing an effect at this level would 
likely require the sensitivity of a future experiment, as 
one might envision at an EIC. 

In the same work [3] it is argued that one could also 
explain the large measured top quark forward-backward 
asymmetry A f FB at the Tevatron [24], if the Z' has a 
fairly large non-diagonal coupling between up and top 
quarks. These particles would certainly generate addi- 
tional contributions to the C\ u and C 2u coefficients, but 
as we have observed this contributions have to examined 
on a model by model basis, since their new Z' bosons are 
not much heavier than the top quark. 

To summarize, we have shown how the future measure- 
ment of Ci q and C 2q will be sensitive to several couplings 
of a leptophobic Z', and will provide interesting informa- 
tion about the flavor structure of these BSM models. As 
an example, we see that a disagreement of future mea- 
surements with the C\ u SM prediction could only be due 
by a Z'tu coupling, whereas a disagreement in C\d could 
hardly be explained by leptophobic Z' . On the other 
hand, a disagreement only in C 2u could probe some of the 
leptophobic Z' models proposed to explain the Tevatron 
W+dijet anomaly. It is important to emphasize that the 
corrections to the C 2q coefficients due to flavor diagonal 
couplings vanish if the Z' does not have axial-vector cou- 
plings to the qq pair. For this reason, models with only 
vector couplings, like the one studied in Ref. [1], cannot 
be probed with PV DIS electron-deuteron scattering. 

More generally, the leptophobic Z' models studied in 
this paper, wherein the coefficients C\ q might not re- 
ceive any significant correction, show nicely how future 
C 2q determinations via Ap 1 ^ are an interesting BSM 
probe, complementary to the precise C\ q determinations 
through atomic PV and PV elastic ep scattering. 
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